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Abstract 

We study unitary random matrix ensembles in the critical case 
where the limiting mean eigenvalue density vanishes quadratically at 
an interior point of the support. We establish universality of the limits 
of the eigenvalue correlation kernel at such a critical point in a double 
scaling limit. The limiting kernels are constructed out of functions 
associated with the second Painleve equation. This extends a result 
of Bleher and Its for the special case of a critical quartic potential. 

The two main tools we use are equilibrium measures and Riemann- 
Hilbert problems. In our treatment of equilibrium measures we allow 
a negative density near the critical point, which enables us to treat all 
cases simultaneously. The asymptotic analysis of the Riemann-Hilbert 
problem is done with the Deift/Zhou steepest descent analysis. For the 
construction of a local parametrix at the critical point we introduce 
a modification of the approach of Baik, Deift, and Johansson so that 
we are able to satisfy the required jump properties exactly. 

1 Introduction 

We consider the unitary random matrix model 

Z-^j^exp{-NTTV{M))dM (1.1) 

defined on Hermitian nxn matrices M in a critical regime where the limiting 
mean density of eigenvalues vanishes at an interior point. It is a basic fact 



1 



of random matrix theory ^3 01] that the eigenvalues of the random matrix 
ensemble (jl.ip follow a determinantal point process with correlation kernel 



n-l 

N-,r/„.\ AT, 



ir„,^(a;, y) = e-T^(-)e"V'fe) Y,PkA^)PkAy)^ (1-2) 

fc=0 

where pk,N denotes the fcth degree orthonormal polynomial with respect to 
the weight e~^^^^^ on M. 

We assume in this paper that the confining potential : M ^ M in 
is real analytic and that it satisfies the growth condition 

' +00 as \x\ — > +00. (1.3) 



log(x^ + 1) 

These assumptions ensure that the mean eigenvalue density -Kn,N{x, x) has 
a limit as ra, oo, n/N — > 1, see e.g. [T7], which we denote by ipvix). 

It is known that tpv is the density of the measure fiy which minimizes the 
weighted energy 

Iv{f^) = jj log j-^^dfx{s)dfx{t) + jv{t)dfx{t) (1.4) 

among all probability measure on M. The measure nv is called the equilib- 
rium measure in the external field V. The fact that V is real analytic ensures 
that the support Sy = supp(/iv/) consists of a finite union of intervals [TH]. 

It is a remarkable fact that local scaling limits of the kernel ()1.2p depend 
only on the nature of the density ipv This has been proved rigorously in the 
bulk of the spectrum for a quartic V in and for general real analytic V in 
Indeed, if ipvix*) > 0, then 

hm -^—-K^,n (x* + — X* + -^^) = K'-'^u, v) (1.5) 
exists, and 

K'^^\u,v) = '^^^^P^. (1.6) 
7r(n — v) 

The scaling limits are different at special points of the spectrum. At edge 

points of the spectrum the density ipv typically vanishes like a square root, 
and then it is known that for some constant c > 0, 

hm -4i7ji^„„„ (x* + X* + = K^'nu. v) (1.7) 
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where 

Ai{u) Ai'{v) - Ai'{u) Ai{v) 



u — V 

and Ai is the Airy function. The Airy kernel is related to the Tracy- Widom 
distribution In ()1.7p we have assumed that x* is a right edge point. For 
a left edge point we change u i-^ —u, v i— ^> —v in the left-hand side of ((TTTj). 
Other special points in the spectrum include 

• Edge points where the density vanishes to a higher order. The possible 
edge point behaviors (at a right end point x*) are 

tPy{x) = c{x* - xf'+^l + o{l)) asx^x*+ (1.9) 

where c > and is a non-negative integer. 

• Interior points where the density vanishes. Then 

i;v{x) =c{x-x*y^{l + o{l)) asx^x* (1.10) 
where c > and k is a positive integer. 

In these critical cases it is believed that the local scaling limit at x* of the 
kernel only depends on the order of vanishing of the density at x* |H|. 

The case where ipv vanishes quadratically at an interior point of Sy, 
that is, the case = 1 in p.lO|) . was considered by Bleher and Its lOj for 
the case of a critical quartic potential V{x) = jx^ + |a;^, with g > and 
t = tc = —2^. Then 



M^) = ^3x'^— - x\ for X e [-2g-'/\2g''/% 

so that ipv vanishes quadratically at the origin. Bleher and Its consider the 
double scaling limit where t changes with n and tends to tc as n ^ oo in 
such a way that n^^^it — tc) remains constant. For the quartic potential this 
is equivalent to considering ()1.1|1 where n,N — > oo, n/N — >• 1, such that 

lim n^/^(--l) (1.11) 
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exists. Bleher and Its gave a one-paremeter family K'^^'^^{u,v; s) of limiting 
kernels, depending on s G M, so that for some c > 0, 

where s is proportional to the value of the limit p.llj) . 

The critical kernels are expressed in terms of so-called ?/^-functions asso- 
ciated with the Hastings-McLeod solution of the Painleve II equation |2H]- 
Consider as in JU] the linear differential equations for a 2-vector (or 2x2 
matrix) = \I'(C; s), 

^^ = A^, 1-* = ^^ (1-13) 

dQ OS 



where 



and 



A-AiC-s)-\ Ae + s + 2q' + 2r^ 



B = B{C-s)=[ \ ^ J. (1.15) 



-C -q 



The compatibility condition for ()1.13|) is that q = q{s) satisfies the Painleve 
II equation q" = sq + 2q^ and that r = r(s) = q'{s). We assume that q{s) is 
the Hastings-McLeod solution of Painleve II, which is characterized by the 
asymptotic condition 

q{s) = Ai(s)(l + o(l)) as s ^ +oo. 

The critical kernels op are given by 



'K[U — v) 

where ( ^2^^'^^ I is the special solution to (|1.13p which is real for real C, 
satisfies 
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and has asymptotics on the real hne 

as C ^ ±00. 

$2(C;s) = -sin (^^c' + <j +0(r') 

If we put 

$^ = $l+i$2^ $2 = $l-i<|.2 (1.17) 

then 

K {u,v,s) = — r (1.18) 

2m[u — V) 

and I ^ 1 is a special solution of the differential equations 

-vE^(C,s)-^ 4^^_2zr 4<2 + ,(, + 2g)J*^^'"^ ^^-^^^ 

and 

The equations p.l9p - p.20p for the V^-functions correspond to the ones used 
by Flaschka and Newell j24j and we will also use those in what follows. The 

vector (^^j is the unique solution of p.l9|) with asymptotics 



as ( ^ 00 uniformly in e < arg ( < tt — e for any e > 0. 

Before discussing our results, we like to point out an integral formula for 
the kernel K'^^^^{u,v; s). If we take a derivative of ()1.16|) with respect to s 
and use p.l3|) and p.l5|) . we get after some calculations 

as TV 

Using the Deift/Zhou steepest-descent method for s —00 as done in ^3] , 
one can show that K'^^^^{u, f ; s) — as s ^ —00, so that we get 



-j [^^{u]a)^^{v;a) + ^'^{u;a)^^{v;a)\da. (1.22) 
^ J —00 



Since $^(C; s) and $^(C; s) are real for real (, formula (|1.22p clearly shows 
that K'^"^(u,u; s) > 0, as it should be. 



2 Statement of results 

It is the aim of this paper to show that the kernel K'^"^{u, v; s) is a universal 
limit. Whenever the limiting mean eigenvalue density ipy vanishes quadrat- 
ically at an interior point, the correlation kernel Kn,N has a double scaling 
limit given by p.l2j) . 

In our Theorem 12.11 below, we use the equilibrium measure ujs of a com- 
pact set S* C M. This is the unique probability measure on S that minimizes 
the logarithmic energy 

J(/i) = jj log j-^^dfx{s)dfx{t) (2.1) 

among all Borel probability measures fi on S. If S* is a single interval [a, b], 
then Us has a density ws given by 

^s(a;) = — . \ =, X e (a, b). 
'Ky [b — x)[x — a) 

If is a finite union of disjoint intervals, say S = IJj=i['^i5 ^j] with aj < bj < 
Oj+i. Then us has density 

wsix) = ^= l^^"^^' xe[j{a„b,), (2.2) 

T^^]X=i\ibj-x){x-aj)\ j=i 

where p{x) is a monic polynomial of degree n — 1 with exactly one zero in 
each of the gaps {bj, cij+i), j = 1, . . . , — 1, see e.g. jHZl Lemma 4.4.1]. Note 
that ()2.2|) has an extension to an analytic function in C \ (M \ 5*°), where 
S° = [Jj{aj, bj), which is a fact that we will use in what follows. 
The following is our main result. 

Theorem 2.1 Let V be real analytic on M such that lim , = +oo. 

Let ipv the density of the equilibrium measure in the external field, and let 
X* be an interior point of Sy = supp{ipv) which is such that 

^y{x*)=^'y{x*)=0, ^'^{X*)>0. 
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Let n, N oo such that the limit 

lim n^/^ (—-l]=L 

exists with L G M. Let Kn^N be the correlation kernel fll.2|) for the eigenvalues 
of the random matrix model . Then there exist constants c > and s G M 
such that 

lim --^K^,N X* + X* + = K^"\u, V- s) (2.3) 

uniformly for u,v in compact subsets o/M. 

Explicit formulas for the constants c and s are 



c = — 

8 

and 



(2.4) 



s = L^^wsA^*), (2.5) 
where Wsy is the density of the equilibrium measure of Sy ■ 

As noted before, Bleher and Its fO] proved ()2.H|1 for the case of a critical 
quartic V. See ^I] for a rigorous expansion of the free energy in this critical 
case. 

Remark 2.2 The random matrix model (jl.lll may be generalized to include 
a spectral singularity at the origin 

Z-Jvl det Mp" exp(-iV Tr V{M))dM, a > -1/2. (2.6) 

If ipv{^) > and n = N ^ oo, then the scaled limit of the correlation kernels 
is a Bessel kernel which involves Bessel functions of order a ± | , see j21 EB] ■ 
In the multicritical case where ipv vanishes quadratically at 0, an analogue of 
Theorem 2.1 is valid. In work in progress we are considering the double 
scaling limit of ()2.6|) and we show that the limiting kernels are expressed 
in terms of the -j/^-f unctions associated with a special solution of the general 
Painleve II equation 

q" = sq + 2q^ - a. 

The main ingredients in the proof of Theorem 12.11 are equilibrium mea- 
sures and Riemann-Hilbert problems. We give some comments on both. 
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-2 2 -2 2 -2 2 



Figure 1: The density of /if for V{x) = ^ — with t equal to 0.8, 1 and 1.3, 
respectively. 



Equilibrium measures. Recall that the equilibrium measure in external 
field V minimizes ()1.4|1 . We need to know how the equilibrium measure /j^v 
in the external field changes as a result of a change in V. The particular 
modification we consider here is 

Vt = jV, t> 0, 

so that Vi = V. Let us put 

Then it is known that tfit and St are increasing as a function of t, see e.g. 
[T^ EHl EH] ■ We also have the Buyarov-Rakhmanov formula |l4j 

1 /•* 

^t = - I t^supp(M.)C?T, (2.7) 

^ Jo 

which expresses the equilibrium measure in the external field as an average 
of equilibrium measures of sets. A consequence of ()2.7|) is that 

= ^Sv, (2.8) 

t=i 

which partly explains why the equilibrium measure Usy plays a role in the 
formula ()2.5j) for s. 

For the case of interest in this paper we have that ipv vanishes at x = x*. 
Then for t > 1, there is a positive density at x*, while for t < 1, x* is out 
of the support of fif For t slightly less than 1, there is a gap in supp(/it), 
see Figure [TJ An asymptotic analysis based on the equilibrium measure fit 
would require a discussion of the two different situations t > 1 and t < 1, as 
is done in 0. 
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Figure 2: The density of i/t for V{x) = ^ — with t equal to 0.8, 1 and 1.3, 
respectively (compare with Figure Q). 



Therefore we found it convenient to introduce a modification of the equi- 
librium problem in external field, which will enable us to treat both cases 
simultaneously. The modification we make is that we do not require the mea- 
sure to be non-negative in a neighborhood of the point x*. For a sufficiently 
small 5o > 0, we consider the problem to minimize 

Iv^i^y) = jj \ogj^^^du{x)du{y) + jj V{x)du{x) (2.9) 

among all signed measures z/ = z/+ — z/~ on M, where v"^ are nonnegative 
measures, such that 

J du = 1, and supp(z/^) C [x* — ^o, a;* -|- (5o]. (2-10) 

We denote the minimizer by ut and we let 5*^ = supp(z/t). 

Then z/j = fit for t > 1, but for t < 1 there is a clear distinction between 
Ut and /it, see Figure 2. However we still have the analogue of ()2.8|) (as we 
prove) 

j^itut) ^^^=ujs,, (2.11) 

What's more, we also have (|2.1ip at the level of densities, that is, if ipt denotes 
the density of Uf, and if x is an interior point of Sy (in particular if x = x*), 
then 

j^itM^)) =wsAx), (2.12) 

where ws^ is the density of the equilibrium measure of Sy, which is what 
we need for the proof of Theorem 12.11 Note that we do not have ()2.12|) for 
X = X* if ipt is the density of fit- 
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Riemann-Hilbert problem. The second main tool for the proof of The- 
orem 12.11 is the characterization of orthonormal polynomials by means of a 
Riemann-Hilbert problem, due to Fokas, Its, and Kitaev [23], and the sub- 
sequent asymptotic analysis of the Riemann-Hilbert problem by means of 
the Deift/Zhou steepest descent analysis of Riemann-Hilbert problems, in- 
troduced in ^JI. , and further developed in QSl I2D1 ■, and other more recent 
papers. The Deift/Zhou steepest descent analysis of the Riemann-Hilbert 
problem consists of a sequence of explicit transformations, which result in 
a Riemann-Hilbert problem that is explicitly solvable in terms of Neumann 
series. (In fact, in this paper, we only need the first term of this series.) The 
critical point x* needs special attention. 

Of particular interest for us is the paper ^ by Baik, Deift, and Johansson 
on the length of the longest increasing subsequence of a random permutation 
of {1, 2, . . . , n}. These authors show that the fluctuations of this random 
variable are distributed according to the Tracy- Widom distribution jSH] in 
the limit as n ^ oo. One of the technical tools in this important paper is the 
asymptotic analysis of a Riemann-Hilbert problem on the unit circle which 
is related to an equilibrium measure (also on the unit circle) whose density 
vanishes at the point —1. This situation is comparable to ours. The authors 
of |S| , see also subsequent papers IH El El E| , construct a local parametrix 
near —1 with the aid of the T/^-f unctions associated with the Hastings-McLeod 
solution of Painleve II. These ^/'-functions satisfy a model Riemann-Hilbert 
problem and the local parametrix is constructed by appropriately mapping 
the model Riemann-Hilbert problem onto a neighborhood of —1 so that it 
satisfies certain desired jump properties approximately. 

We follow the approach of P] but we introduce a modification in the con- 
struction of the local parametrix so that it has the desired jump properties 
exactly, in contrast to |B] where the desired jump properties only hold ap- 
proximately. The fact that we have the exact jump properties simplifies the 
arguments considerably and we feel that this is also a main contribution of 
the present paper. 

Outline of the rest of the paper. In Section 3 we collect the necessary 
facts about equilibrium measures. In particular we study the modified equi- 
librium problem with external field in some detail. In Section 4 we discuss 
the Riemann-Hilbert problem satisfied by the V'-functions. Here we follow 
|24j . Then in Section 5 we state the Riemann-Hilbert problem for orthogonal 
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polynomials, discuss the relation with the correlation kernel Kn^Ni and per- 
form the transformations in the steepest descent analysis. Finally in Section 
6 we give the proof of Theorem 12.11 

3 Equilibrium Measures 

As explained in the previous section, we consider a modification of the equi- 
librium problem where we drop the non-negativity condition in a small neigh- 
borhood of X* . We take 5o > sufficiently small so that 

V'y(x) > for all xe[x* - Sq, x* + 5o] \ {x*}. (3.1) 

and we use ut to denote the signed measure that minimizes ()2.9p under the 
conditions (j2.1U|) . We define 

St = snpp{ut). (3.2) 

The existence and uniqueness of ut follows as in |36j . 
Let 

f/^(x)= / log^^dz/(y) 
J \x-y\ 

be the logarithmic potential of u. Then standard arguments of potential 
theory [121 ES! show that Ut is the unique signed measure satisfying ()2.9j] 
with the property that 

2U'''{x) + jV{x) = it, X esnpp{ut)U[x* -6q,x* + 6o], (3.3) 

2f/^*(s) + iy(x) > £t, xeR. (3.4) 

for some constant it- 

As before we use /if to denote the equilibrium measure in the external 
field Vt. This is a probability measure that satisfies for some constant it, 

2f/^*(x) + jV{x) =it, xe supp(/iO, (3.5) 

2U^''{x) + jV{x) >it, xeR. (3.6) 

It is known that tfit and supp(/i() are increasing with t > 0, see |14| I36j. For 
t > 1, we have supp(/iy) C supp(/it) which implies that ut = fit in view of 
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(|3.1|) and the variational conditions (|3.3|) - (j3.6p . For t < 1, we have that x* 
is outside the support of /it, and in fact the strict inequahty in (|3.6p holds 
for X = X* ii t < 1. 

In general there is the following inequality between fit and ut. 

Lemma 3.1 For every t > we have 

/it < lyt- (3.7) 

Proof. Let X = Ut - /i*. From (jSISl), (D, (ESI), and dHIHl) it follows that 

2U^{x) <et- it, for X G St, (3.8) 
2U^{x) >it- it, for X e supp(/it). (3.9) 

The potential is subharmonic on C \ supp(A^) and since J dX = 0, it is 
subharmonic at infinity as well. By the maximum principle for subharmonic 
functions I36j. the maximum of is attained in supp(A^) only. Since 
supp(A+) C St we then have by ()3.8p that equality in ()3.9p holds for every 
X e supp(/it), and so supp(/it) C supp(A"'"). This implies (j3.7p . □ 

It follows from ()3.7|) that supp(z/j^)nsupp(/ij) is empty. Since for t slightly 
less than 1 a gap opens in supp(/if), which depends continuously on t, see 
jS21, it follows that for any given 6 G (0, 6o), there is to < 1 such that 

supp{h'^) C [x* — 6,x* + 6], for t > to- 

This shows that for t < 1 sufficiently close to 1, the definition of ut is inde- 
pendent of the choice of 6o. 

A very useful fact is that for real analytic V, say V is analytic in a 
neighborhood V of the real line, the measures fit have densities ipt which 
can be expressed in terms of the negative part of an analytic function in V. 
Indeed, if 

then it was shown in jTH] that 

Mx) = -\/^), (3.11) 

where = max(0, —qt) denotes the negative part of qt- A consequence of 
(I3.11|) is that supp(/it) is the closure of the set where qt is negative. 
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The arguments of [TH] can be readily extended to the signed measures Ut, 



provided that supp(z/j ) C [x* 



2t 



6, X* + 6] for some 6 < 60. We define 
V'{z) - V'{y) 



z-y 



zeV. 



(3.12) 



and then the following holds. 

Proposition 3.2 There exists to G (0, 1) 
signed measure Vt has a density tpt o^nd 

qt{x) = - [vr^t(x)]^ , 



such that for every t > to, the 



for X G Sf. 



In addition we have 



Proof. This follows as in Proposition 2.51]. 



forx eR\St 



(3.13) 

(3.14) 
□ 



Obviously, for t > 1 we have qt = qt- Since ipi = ipy has a double 
zero at x*, we see from (I3.13p that qi has a zero at x* of order four. For t 
slightly bigger than 1, this fourth order zero splits into two double zeros in 
the complex plane away from the real axis. 

For t slightly less than 1, there is a difference in the behavior of the zeros 
of qt and qt- Indeed, qt has two simple real zeros near x* which are endpoints 
of the support of fit, and in addition there is a double real zero in between 
them. On the other hand, we have that qt has two double real zeros near x*, 
which are endpoints of the support of . 

The fact that —qt has only double zeros near x* allows us to take an 
analytic square root, and in view of ()3.13p we choose it so that for z & St, 



Mz) = -i-qtiz)Y/\ 

TT 



(3.15) 



where the sign of the square root at x = x* is taken negative if t < 1 and 
positive if t > 1. The right-hand side of ()3.15|) has an analytic extension to 
a neighborhood of x*, which is independent of the value of t > to- Thus ipt 
has an analytic extension to a fixed neighborhood of x*, which will also be 
denoted by ipt- 

Recall from the discussion before the statement of Theorem 2.1, that the 
density ws^ of the equilibrium measure of Sy also has an analytic extension 
to a neighborhood of x*, which we also denote by wsy- The remaining part 
of this section is devoted to the proof of the following proposition. 
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Proposition 3.3 We have 

Jim -j—^ = wsy{z) (3.16) 

uniformly for z in a neighborhood of x* . 

We start with a lemma which contains a weaker form of ()3.1(ij) . 
Lemma 3.4 We have 

uosy (3.17) 



dt 

where uJsy is the equilibrium measure of Sy- 



Proof. Buyarov and Rakhmanov J3] proved (for a very general class of V) 
that 

and 

hm ^^^i^ = c^,, (3.19) 

where us^, is the equilibrium measure of the set 

5^ = {a; G M I 2f/^'^(x) + V{x) = £i}, 

see fl3.5p with t = 1. For real analytic V, the sets Sy and Sy differ by an at 
most finite number of points, so that = ^s^,- Thus by fl3.18p and ()3.19p 
we have 

Now write 

In view of fl3.2()|l it suffices to prove that 

||z/i — jitW = o(t — 1) as t — 1, (3.21) 

where the norm denotes the total variation of a signed measure. We may 
assume that t < 1. Because of ()3.7|1 we have 



- vt II = y A'^t y = 2 y dv-. (3.22) 
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For if: < 1, the support of z/^ is contained in an interval [x* — 5{t), x* + S{t)], 
with 

=0((l-t)^/2) ast-^1-, (3.23) 

see 1221 Lemma 8.1(iii)]. From (j3.15|) and the fact that qt has two double 
zeros in [x* — 6(t),x* + 6(t)] which are the end-points of the support of , 
we then easily get that 

= 0{5{ty) = 0(1 - 1) as t ^ 1 - . (3.24) 

Combining (IT^ and (jT^ . we find / du^ = O ((1 -tf/^) as t 1-, 
which by ()3.22j) implies ()3.21|) . This completes the proof of the lemma. □ 

We now give a characterization of ipv, which will be of use in the proof 
of Proposition 13.31 

Lemma 3.5 For x E Sy we have 

M^) = / '^^^^^^i^s.. (.), (3.25) 

2Ti^wsy[x)J x-y 

where ws^ is the density of usy ■ 

Proof. By the Sokhotski-Plemelj formulas, see e.g. we have that 
F{z) = I ^i^^ds, foTzeC\ Sv, 

JSv 

satisfies 

F+(x)+F_(x) = V\x), xeSv, (3.26) 

F+(x)-F_(x) = -2milJv{x), x e Sy (3.27) 

Let Sv = Uj'=i[%5^i] set 

/ n \ 1/2 

R{z) = [jliz - h){z - a,) j , zeC\Sv, 

where the square root is positive for z > bn- Using ()3.26p and the fact that 
R+{x) = —R-{x) as a; G Sy, we see that 

p{x)F+{x) _ p{x)F_{x) ^ p{x)V'{x) 
R+{x) R^x) R+{x) ^ ^ 
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for any polynomial p. Suppose p has degree at most n — 1. Then (j3.28|) and 

as z — * cxD imply that 



the fact that ^^^j^ 

R{z) 

p{z)F{z) 



R{z) 



1 

27ri 



p{s)V'{s) 



Sv 



RAs) 



-ds. 



which we can rewrite for z E V as 

p{z)F{z) _ 1 

1 

~ 2Tii 



R{z) 



V'is) - V'{z) p{s) ^ V'{z] 



zeC\Sv 
p{s) 



Sv 



RJs) 



2m 



Sv 



R^(s) s — z 



-ds 



V'{z) - V\s) p{s) ^ V'{z) p{z) 



Sv 



RJs) 



2 R{z) 



(3.29) 



by contour integration. Now we use ()3.29j) and ()3.27|1 to obtain 



F+{x) - F_{x) 
-27ri 



1 R+{x) 
p{x) 



V'{x) - V'{y) p{y) 

Sv ^-y R+iv) 



dy, 



and this holds for any polynomial p of degree at most n — 1. Since we know 
by (I2I2I) that 

i p{x) 



u^SvK^) = — 



X G 5'v 



TT R+{x) ' 

for some monic polynomial p of degree n — 1, we get ()3.25|) . 

Now we can give the proof of Proposition 13.31 
Proof of Proposition 13.31 By fj3.12|) we have for z G V, 



□ 



V'{z) 



V'iz) - r{y) 



z-y 



d{tut){y) 



so that in view of (|3.17|) we find, 



V'{z) - V'{y) 



z-y 



d^Sviy)- 



Thus by fl3.15|) for 2; 7^ x* in a neighborhood of a;*, 



t=i 



dt TT 



1 

2^ 
1 

27r2 ^y{z) 



1 



t=i 
Viz] 



V'{y) 



z — y 
V'{z) - V'{y) 



duJsviy) 



z-y 



duJsviv)- 



(3.30) 
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By (|3.25|) . the right-hand side of (j3.3U|) is ujsv{z) in case z ^ x* is real, and by 
analytic continuation it continues to be {z) in a punctured neighborhood 
of X*. Thus ()3.16|) holds for z ^ x* m. a, neighborhood of x*, and then it 
easily follows for z = x* as well. □ 



4 Riemann-Hilbert Problem Associated with 
the Painleve II Equation 

In this section we recall the RH problem associated with the second Painleve 
equation, see [211 EOl 123 13 121] • We transform this RH problem into a RH 
problem that will be used in the next section to construct a parametrix 
around x* . 



4.1 ?/;-functions Associated with Painleve II 

Consider the differential equation of ()1.19|) . 

dC^^' ~ V 4Cg - 2ir AiC + i{s + 2g2) ) ^^"^^ 

where \l/ is a 2 x 2 complex matrix- valued function and s, q and r are con- 
sidered as parameters. All solutions of 1)4.11) are entire functions of C,. 
For j = 1, . . . , 6, let Sj be the sector 

= {C G C I ^^TT < argC < ^^vr}. (4.2) 



There exists a unique solution \E'j of equation ()4.1)1 so that 

*j-(C)e'^^^'+'^^"=' = / + 0(C"') (4.3) 

as C ^ cxD in the sector Sj. Here we use ~ (^0 '^l) denote the third 

Pauli matrix. There exist complex values a^, j = 1,...,6 (called Stokes 
multipliers) so that 

vI/,+i(C) = vl/,(C)A„ forj = l,...,5, vl>^(C) = vl>g(C)A6. (4.4) 

with 

A,= Q °) if J is odd, (4.5) 
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and 

Aj = if j is even. (4.6) 

Furthermore, we have 

= ttj and 0102^3 + fli + 02 + % = 0. (4-7) 
Now define the rays 

r, = {CI arg C = ^^^} for J = 1, ... , 6, (4.8) 
oriented away from the origin, and the matrix-valued function 

^{0 := ^,{0, for C e S,. (4.9) 
Then \1/ satisfies the following RH problem 

(a) is analytic on C \ IJ^=i ^j- 

(b) ^+ = ^^Aj on Tj, 

(c) ^(C)e^(t^'+<)'^» = / + 0(C-^) as C ^ 00, 

(d) \1/ is bounded near 0. 

The Stokes multipliers aj depend on s, q and r. An isomonodromy de- 
formation is a variation of these parameters such that the Stokes multipliers 
remain constant. Flaschka and Newell [23] showed that the isomonodromy 
deformations are given by the Painleve II equation q"{s) = sq{s) + 2q'^{s) 
and r(s) = q'{s). 

Any solution of the Painleve II equation is a meromorphic function with 
an infinite number of poles. We write \1/(C; s) for the \I' function ()4.9|) with 
parameters s, q = q{s), r = r{s), where q is the Hastings-McLeod solution 
and r(s) = q'{s). Let V be the set of poles of q. Then \I'(C; s) is defined and 
analytic for C G C \ IJj=i fo^ s E C\V. It is known that there are 

no poles on the real line [2H]- The Stokes multipliers corresponding to the 
Hastings-McLeod solution are 

ai = 1, a2 = 0, 03 = —1. 
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Figure 3: Jumps for s). 

Thus \1/(C; s) is analytic across the imaginary axis. We reverse the orientation 
of Fs and F4, and we define 

El = Fi U F3, and E2 = F4 U Fg. 

Then s) solves the following RH problem, see also Figure 0] 

(a) \1/ is analytic on C \ (Si U S2), 

(bl) ^+ = ^ 
(b2) ^+ = ^ 

(c) ^(C; s)e*(l«'+<)'^=' = / + OiC^) as C 00, 

(d) \1/ is bounded near 0. 

The RH problem has a solution if and only if s G C \ P. The properties (c) 
and (d) are valid uniformly for s in compact subsets of C \ P. 
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1 
1 1 



1 -1 
1 



on Si, 



on S 



2, 



(4.10) 




Figure 4: Jumps for s). 

4.2 Transformation of the RH problem 

We modify the RH problem for \E' so that it resembles the RH problem that 
we will need locally near x*. 

We introduce an additional parameter G M and define 



gi0<T3^(^. 3)gi(3C^+<)<x3g-ie<73 for ImC > 0, 

-1 

1 



MiC,s,9) { ^ieas^^^.^y^i^s+sCWs^-ieaJ^ 1| for ImC <0. 



(4.11) 

For any given s E C \V and ^ G M, we then have that M is defined for 
C G C \ (M U El U S2), see Figure and satisfies the following RH problem 

(a) M(C; s, 9) is analytic for C G C \ (M U Si U S2) and s G C \ P. 

(bi) M^ic,s,e) = M4c,s,e) (^^,,^,J_,,^_,) °) for c e Si, 

(b2) M+(C; s, 9) = M„(C; s, 9) (^_°^ for C E R, 

(b3) M+(C; s, 9) = M_(C; s, 9) (^^_2,(|c^3+,^_,) °) for C G S^. 

(cl) M(C; s, 0) = / + 0(C~^) as C ^ 00 in the upper half plane, (4-12) 
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Figure 5: Jumps for M{(; s, 6). 




as C 



oo in the lower half plane, 



(4.13) 



(d) M((^; s, 6) is bounded for ^ near 0. 

The properties (cl), (c2), and (d) hold uniformly for s in compact subsets of 
C \ P and for G M. 

5 Riemann-Hilbert Analysis 

The proof of Theorem 12.11 is based on the steepest descent analysis of the 
Riemann-Hilbert problem for orthogonal polynomials. 

Since the main point of the present discussion is the treatment of the 
critical point x*, we will restrict ourselves to the one-interval case. We also 
assume that there are no other singular points besides x* . Thus Sy = [a, 6] 
and ipv vanishes like a square root at a and 6, and quadratically at x* but 
at no other points of Sy. In addition, we assume that the inequality in the 
variational condition (|3.6|) is strict for x G M \ [a, 6] . We then have that 



for certain at < ht for every t close to t = 1. Note that at is an increasing 
and ht a decreasing function of t. 

We will comment below on the modifications that have to be made in the 
multi-interval case, see Remark 16.11 



St = supp(z/t) = [at,bt] 



(5.1) 
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5.1 RH problem for Orthogonal Polynomials 



For each n and A^, we consider a Riemann-Hilbert problem, introduced by 
Fokas, Its and Kitaev |2Sll2ni- We will look for a 2 x 2 matrix- valued function 
Y = Yn,N (we drop the subscripts for simplicity) that satisfies the following 
conditions: 



(a) Y is analytic in C \ M, 

(b) Y+{x) = r_(x) 

(c) Y{z) = {I + 0{z~')) 



I ^-NV{x) 
1 



for X G 



(5.2) 



2" 
Z-' 



as z — >• oo. 



Here V+(x) (resp. Y^{x)) denotes the limit as we approach x G M from the 
upper (resp. lower) half-plane. The RH problem possesses a unique solution 
which is given by 



Yiz) 



( 



\ 



-27rm„„i^Arp„_i,Ar(2) -K„_i,Ar 



-NV{s) 



-ds 



\ 



s — z 



Pn-l,N[S)e' 



-NV(s) 



-ds 



(5.3) 



s — z 



where Pn,N{x) = Kn,NX^ + ■ ■ ■ denotes the orthonormal polynomial as before, 
see also pTf l!TT]. 

The correlation kernel ()1.2|) can be expressed directly in terms of the 
solution of the RH problem. Indeed, by the Christoffel-Darboux formula for 
orthogonal polynomials 

f ^_ -i^V(x) -^V{y)'^ri-l,NPn,N{x)Pn-l,N{y) -Vn^l,N{x)Pn,N{y) 
-l^n,N{X,y) —62 e 2 

Kn,N X — y 

which involves the orthogonal polynomials of degrees n and n — 1 only. Using 
()5.3p and the fact that det Y{z) = 1 for every 2; G C \ M, we then get 



i^„.^(x,y)=e-f^(^)e-f^(^) 



2TTi{x — y) 



(0 l)Y-\y)Y4x)(^^. (5.4) 



Applying the Deift-Zhou steepest-descent method to the RH problem will 
allow us to find the asymptotics of Kn^N as given in Theorem 12.11 
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We note that an expression like (j5.4j) was first given in the context of 
random matrices with external source [U ^1 where the eigenvalue corre- 
lation kernel is given in terms of the solution of a 3 x 3 matrix valued RH 
problem, and the expression similar to ()5.4|) was found to be very convenient 
for asymptotic analysis. Also in the present 2 x 2-case we find it helpful to 
work with (|5.4|) . 



5.2 Normalization of the RH problem 

In the first transformation we normalize the RH problem at infinity. A 
standard approach would be to use the equilibrium measure fit in external 
field Vt where t = n/N. If we would do this for the case where t < 0, we 
would have an equilibrium measure with a gap in the support around x*, and 
an annoying consequence is that the equality in the variational conditions is 
not valid near x*. For this reason we have introduced the signed measures 
in Section 2 and we will use these measures now to normalize the RH 
problem. 

We let t = n/N and assume t is sufficiently close to 1 so that (j5.1|) holds. 
As before we use 

P (5.5) 

ax 

to denote the density of Uf. Then we define the (^-function 

gt{z) = / \og{z - s)dvt{s) = / \og{z - s)ilJt{s)ds, (5.6) 



where we take the branch cut of the logarithm along the negative real axis. 
Then the following properties of gt are easy to check. 



i) e^*^^^ is analytic in C \ [at, b, 



=9t(z) 



z + oa 



n e^'^^ = z + (J - as 2; — > cxD. 



(iii) gt+{x) - gt_{x) = 2Txi i)t{s)ds for a; G M, 

(iv) gt+{x) + gt_{x) - \V{x) + < for x G M \ [a*, h], 

(v) gt^{x) + gt^{x) - \V{x) + £t = for a; G [at, h]. 
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In terms of the analytic function g^, see (|3.12p and Proposition 13.21 we 
have 



rx 

I 1/2 

gt+{x) - gt_{x) = -2 / {qt{s))i ds, for x e [at,bt], (5.7) 

Jbt 



and 



gt+ix) + gt_{x) - \v{x) +et = -2 [ {qt{s)y^^ds for x > (5.8) 

Jbt 

gt+{x) + gt^{x) - \v{x) + it = -2 [ {qt{s)y/^ds for x < at (5.9) 

^ J at 

where {qtis))^^"^ is analytic for s G V \ [at,bt] and that square root is taken 
which is positive for large real s. We define 

Vt{z) = f\qt{s)f"ds (5.10) 
Jbt 

and 

Uz) = !\qt{s)fl'ds (5.11) 



which are defined and analytic in the neighborhood V of the real line where 
V is analytic with cuts along (— oo, ht) and (oi, +cxd), respectively. Note that 
qt has simple zeros in at and ht and only double zeros in V \ [at, ht], so that 

1/2 

is indeed analytic there. It is possible that qt has double real zeros, and 
that qt has sign changes in {—oo, at) or {ht, oo). However, by ()5.8p and ()5.9p 
we have that v't(x) > for x > ht and 'ft{x) > for x < aj, since we are in 
a situation with strict inequality in ()3.4|) for x > ht and x < at- 

We will now perform the first transformation of the RH problem: we take 
t = n/N and define 

T{z) = e^^*"3 y(2)e-"9tW-3g-t^t<x3 for ;z e C \ R. (5.12) 

Using the jump condition ()5.2|) of Y and ()5.12p . we easily check that T^{x) = 
T^{x)Jt{x) for X G M, where 

/i „-NV{x)\ 



g-n(9t + (a::)-gt_(2.')) gn(3t + (a;)+gt _ (x)- i V'(a;)+^t) 
Qnigt + ix)-gt-ix)) _ 



(5.13) 
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Because of the properties (j5.7|) - (j5.9|) and the definitions (|5.1(J|) - (j5.11|) . we see 
that the jump matrix Jt has the following forms on the respective intervals 
(at,6j), (fei,oo), and (-00, a^), 

Jr = I Q ^1 on (6t,oo) 
Jr = ( Q ^1 on (-00, at). 

Thus T is the unique solution of the RH problem 

(a) T is analytic in C \ M, 

(b) r+(x) = T_(x)Jt(x) for x G M, 

(c) r(^) = / + 0{z-^) as ^ ^ 00. 

5.3 Opening of the lens 

The jump matrix on the interval (a^, ht) has the factorization 




(5.14) 



Figure 6: The contour after the opening of the lens. 

Now we open lenses around the intervals (a^, x*) and (a;*, ht) as shown in 
Figure El Let Ci be the upper lips of the lenses and C2 the lower lips, with 
orientation as in Figure IHl We open the lenses in such a way that they are 
fully contained in V, the region of analyticity of V and qt. In addition, we 
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can take Ci and C2 in such a way that Rey^t < on Ci and C2, with the 
exception of a neighborhood of x* if t < 1. This effect is due to the fact that 
Vt has a negative density near x* ii t < 1. However, if t < 1 increases to 
one, the exceptional neighborhood shrinks to a point. It follows that for any 
given 5 > 0, there is a constant 7 > such that for t sufficiently close to 1 
we have 

Rev9t(z) < -7 < 



for all z G Ci U C2 with min(|z 
We define 

T 
T 



X 



\z — a\,\z — h\) > 5. 



S 



T 



outside the lenses, 

in upper parts of the lenses, 

in lower parts of the lenses. 



(5.15) 



The RH problem for T and the factorization ()5.14|) imply that S solves the 
following RH problem: 



(a) S is analytic in C \ (M U Ci U C2) , 

(b) S+{z) = S4z)Js{z) for 2 G M U Ci U C2, 

(c) S{z) = I + 0{z^^) as z 00, 

(d) S remains bounded near at, bt, and x*, 
where Js is given by 



(5.16) 




(5.17) 



on {at, bt), 
on C1UC2, 
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The jump matrices on R \ (oj, bt) and on Ci U C2 tend to the identity 
matrix as n — > 00 and t 1. Ignoring these jumps, we find the parametrix 
for the outside region. Uniform convergence breaks down in neighborhoods 
of at, bt, and x*, so that we will also need to construct local parametrices 
near those points. 

5.4 Parametrix away from special points 

The outside parametrix S°° = solves the following RH problem 

(a) S°° is analytic in C \ [at, bt], 

(b) S^ = S^(^^^ on {at,bt), 

(c) S^{z) = 1 + 0{z'^) as z ^ 00. 
As in 1201 it has the solution 



^r(^) = ( /E£)J(£)li ) ' ^eC\[a„6,], (5.18) 



2i 



where 



1/4 

I3{z)=l3t{z)=(^-^] , zeC\[at,bt]. (5.19) 



z-bt 



z-at 



Note that S'^ depends on t. 



5.5 Parametrix at edge points 

At edge points at and bt the density ipt vanishes like a square root. This 
allows the construction of local parametrices near at and bt with the use of 
Airy functions. We will not give the details, see jTTl UHl 1201 • 

5.6 Parametrix at critical point 

In a neighborhood A = {z G C | |a;* — 2;| < 5} oi x* , we want to construct a 
parametrix P with the following properties 

(a) P is analytic in A \ (M U Ci U C2), 
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(b) P+ = P^Jsiz) on (R U Ci U C2) n A, 

(c) P{z) = (/ + 0{n-^/^))S^{z) as n ^ 00, t ^ 1, uniformly for z G ^A, 

(d) P{z) remains bounded for z near x*. 

We seek a parametrix P near x* in the form 

P{z) = Et{z)M (n^/'fiz); n''^htiz),net) (5.20) 

where Et is analytic in A, / is a conformal map from A to a neighborhood 
of 0, St is analytic in A, and Of is a real constant. Recall that M is given by 
(gUI). 

In view of the jump properties of M, we seek /, st, and 9t so that for 

ze A, 

-f{zr + s,{z)fiz)-e,]=i iflm.>0, 

Since ipt+ = — ft-, the right-hand side of ()5.2ip does indeed define an analytic 
function in A. 

The conformal map / will not depend on t. Recall that qy has a zero of 
order four at x*, and that it is negative on (a, b)\{x*}. Thus f^, {—<iv{y)Y^'^dy 
is analytic in A and has a zero of order three at x* . We may assume there 
are no other zeros in A. Then we can take a third root and define 



-qv{y)f"dy 



1/3 



2 e A. (5.22) 



Then / is analytic in A with 



f{z)=r-^i^\'\z-x*) + 0{{z-xn-) 



= c^''\z-x*) + 0{{z-x* f) , asz^x* (5.23) 

where c is given by ()2.4j) . Taking smaller A if necessary, we then have that 
/ is indeed a conformal map on A. Note that f{x*) = and that / is real 
and positive on {x*,x* + 6) and real and negative on (x* — 6,x*). We still 
have some freedom in opening the lenses. We take Ci and C2 in A so that / 
maps them to the rays where M has its jumps. That is, Ci fl A is mapped 
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into arg^ = 7r/6 and arg^ = Svr/G, and C2 fl A is mapped into aYg( = —n/G 
and arg(^ = —5n/6. 

Having / with f{x*) = 0, we take z = x* in ()5.2H) and we see that we 
should take 

et=tipt4x*) = -iipt_ix*). (5.24) 
Then 6t is real and it is also given by 



et = - r i-qt{y)Y^'dy 
Jbt 



and so 



etTm{z)= I {-qt{y)f'^dy, for±lmz>0. (5.25) 

Having / and 6t we finally take St{z) so that 1)5.211) holds, that is, 

4 



s,{z)f{z) = -^f{zf + etTi^t{z) 



{{-qt{y)Y" - {-qv{y)Y'') dy (5.26) 



X* 



where for the last line we used fl5.22|) and ()5.25p . Since the right-hand side of 
()5.26|) is analytic in A and vanishes for 2; = x* we can divide by f{z) (which 
has a simple zero at 2; = x*) and obtain an analytic function in A. Then 
St{z) is real for real 2;, and Si{z) = 0. 

So the above construction yields analytic functions /, St and a constant 
6t so that ()5.21|) holds. Then we can define the parametrix P by ()5.2U|) . 
provided that n^^'^Stiz) for z G A stays away from the poles of the Hastings- 
McLeod solution of Painleve II, since for s G P we have that M{(; s, 6) is 
not defined. For given values of n and t, we can take A small enough so 
that this is indeed the case. However, we want to let n — >■ 00, t — >■ 1 so that 
n^/^[t — 1) — i> L, and work with a neighborhood A that is independent of n 
and t, although it may depend on L. 

Note that by and TT^ 



st{z)f{z) = vr / (My) - My))dy. (5.27) 

J x* 

Because of f|3.16|) we have 

My) - i^viv) = (t - i){wsAy) - My)) + o{t - 1), 



29 



as t — s> 1, uniformly for y in a neighborhood of x*. Using this in (|5.27|) we 
get that 

st{z)^^ = Ti{t - 1)^- / {ws,-{y) - My))dy + o{t - 1) 



= 7r(t - + {t- l)0{z - X*) + o{t - 1), (5.28) 

where o{t — 1) is uniformly in z as t — ^ 1, and 0{z — x*) is uniformly in t as 
z —>■ X*. By ()5.23p and ()5.28p we then also have 

st{z) = 7rc^/3(t - l)wsyix*) + {t- 1)0(2 - X*) + o(t - 1), (5.29) 

where again o(t — 1) is uniformly in 2; as t — 1, and 0{z — x*) is uniformly 
in t as z — > X*. Then if n^^^it — 1) — > L we get 

n^l^st{z) = nc-^^^Lws^x*) + 0{z - x*) + o(l) 

= s + 0(2-x*) + o(l), (5.30) 

where we used the definition ()2.5|1 of the constant s. Since there are no poles 
on the real line, we can find a neighborhood A of x* such that n'^^^St{z) ^ V 
for all z G A if n is large enough and n^^^(t — 1) — >■ L. Note that A depends 
on L, but not on n and t. Then 

M{n'/^f{zy,n^/ht{z),net) 

is well-defined for 2; in a fixed neighborhood A of x*. 

Finally, we define in such a way that the matching condition at dA is 
satisfied. We do this by defining 

S^{z) for Im^ > 0, 



-1 



so that Et{z) is analytic near x* by the jump property (j5.16|) of 5*^°°. Because 
of the asymptotic behavior of M{(; s, 9) in ()4.12j) and ()4.13|) (which is valid 
uniformly for s away from the poles and for 6* G M) we have the matching 
condition 

P{z) = Sr{z)il + 0(n-V3)) ^ ^ 0{n-'/'))Sriz), 

uniformly for z G dA. This completes the construction of the parametrix P 
in the neighborhood of x*. 
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5.7 Third Transformation 



We set 



R{z) 



S{z)P ^{z) for z in disks around x*, a, and b, 
S{z){S^)^^{z) for z outside the disks. 



(5.32) 



Since 5" and P have the same jumps inside each of the disks, and S and S'^ 
have the same jumps on [at,bt], R has jumps on a contour F as shown in 
Figure H 



Figure 7: The contour F after the third and final transformation. 

R solves the following RH problem: 

(a) R is analytic in C \ F, 

(b) R+{z) = R^{z)Jr{z) as z G F, 

(c) R{z) = I + 0(z"^) as z oo, 

(d) R is bounded, 

for certain jump matrices Jr, which as n ^ oo, t ^ 1 so that — 1) — > L, 




satisfy 




I + 0(n-i/3) 
/ + 0(e-T") 



for z in circles around a and b 



for 2; G 9A, 

for some fixed 7 > elsewhere on F. 



As in |17| IT^ it now follows that 



mz)-i 



00 



0{n 




(5.33) 



uniformly for z G C \ F. 
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6 Proof of Theorem 12.11 

Now we are ready for the proof of Theorem 12 .11 We start with the expression 
()5.4j) for the kernel Kn,N- After the transformation ()5.12|) . we find that for 

x,y e {at,bt), 

i^n,^(a:,y) = ^-^^(0 e-'^^^(y^)T^\y)T4x)[^ ^ j. (6.1) 

Using formula ()5.15|1 for S in the upper parts of the lenses, we get for x,y & 

{at A), 

(6.2) 

Assume that x and y are inside the disk A around x* . Then by ()5.32|1 , ()5.2()|1 , 
dmH), and (IK:^ we get 

S+{x) = i^(x)E^(x)e™^*"3^+(n^/V(a;);n2/3si(x))e"'^'+(")"^ 

and so 

j = i?(x)E,(a;)e^"^'-vl>(^V3;(^). n^^ix)) . (6.3) 

Similarly 

= (-1 l)^-\n'/'f{yy,n'/'st{y))e~'"'^'^-E-\y)R-\y). (6.4) 
Now we fix u and v and take 

^ = ^* + / ^1/3 y = X* + ^ , , (6.5) 

so that for n large enough, x and y are inside the disk around x*, so that 
(Q, O, and (jn3D hold. Then it follows from (PT^ and (jHSD that 

n^^^f{x) — > M, and n^^^fiy) v as n ^ oo. (6.6) 

From (j5.3(jp and (j6.5j) we get 

n^^^St{x) s, and n^^^Stiy) ~^ ^ (6.7) 
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as n — i> oo, t ^ 1 such that n^/^(t — 1) L. Furthermore, from (|5.33|) . (|6.5p . 
and the fact that R is analytic near x*, we get 

R-\y)R{x) = 1 + (i^) =1 + . (6.8) 

From ()5.18|) . ()5.3H) . and ()6.5p we easily get 

Er\y)E,{x) = I + o(^^'^ (6.9) 

as n — s> oo. The constants implied by the 0-symbols in (16. 8|) and (j6.9p are 
independent of u and v, when u and v are restricted to a compact subset of 
M. Combining ()6.8|) and ()6.9|) we get that 

e-'''^"''E-\y)R-\y)R{x)Et{x)e''''''''^ =1 + (^^^^^ , (6.10) 

since 6t is real. 

Then multiplying (j6.4j) and (j6.3j) and letting n ^ oo, t —>■ 1 such that 



n 



(t - 1) L, we get by using (Q, (EH), (Q, (EH), and (jFTTnil that 
lim \ Kn^N (x* + " X* + 



= --i r(-l l)^-\v;s)^iu;s)(]), (6.11) 

uniformly for -u and v in compact subsets of M. 
By dTTTll . dO)) and (jO)) we have that 



In view of the jump fl4.10|) satisfied by \1'(C; s) on Si we have that 

ft(C ^)) ^ (i) ^ ^ ^ ^^-^^^ 



33 



Since det \1/ = 1, we also get (after simple calculation) 

(-1 i)M/-i(C;s) = (-<f2(C;s) $i(C;s)) force S,us,. (6.13) 

Note that 5*1 U 5*4 includes the full real line, so that we can take ^ = m in 
()6.12|) and ( = v in ()6.13p which we use in ()6.11|) to obtain ()2.3|) . This 
completes the proof of Theorem I2.ll 

Remark 6.1 The above proof of Theorem 12. II was given under the assump- 
tion that Sv consists of one interval, and that there are no other singular 
points except for x*. Here singular point refers to the classification of [201 
according to which there are three types of non-regular behavior for a real- 
analytic external field V, see also [S2I- The singular points of type I are 
points in M \ where equality in the variational inequality ()3.(ij) with t = 1 
holds. Singular points of type II are interior points of Sy where the density 
ipv vanishes, and singular points of type III are edge points of Sy where ipv 
vanishes to higher order than a square root. Here we indicate briefly the 
modifications that have to be made if these assumptions are not satisfied. 

If Sv consists of more than one interval, then the main complication is 
that the construction of the outside parametrix S'j^ is more complicated, 
since it uses 9-functions as in [201 Lemma 4.3]. It should be noted that S'j^ 
will also depend on n. As a result it will follow that Et as defined in fl5.3ip 
also depends on n. However this will not effect the asymptotic behavior (|6.9|) 
as n — i> oo, so that the above proof goes through. 

In case there are other singular points, we have to construct special local 
parametrices around each of them. For our purposes, we only need to know 
the existence of those parametrices, which is established in |2m Section 5]. 
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